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Abstract 

The rate for the photon emission accompanying orbital IS electron capture by the atomic nucleus is 
recalculated. While a photon can be emitted by the electron or by the nucleus, the use of the length gauge 
significantly suppresses the nuclear contribution. Our calculations resolve the long standing discrepancy of 
theoretical predictions with experimental data for A J = 2 forbidden transitions. We illustrate the results 
by comparison with the data established experimentally for the first forbidden unique decays of 41 Ca and 
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I. INTRODUCTION 



Orbital electron capture by the atomic nucleus (EC) is one of nuclear decay modes governed 
by weak interactions, a common and well known type of radioactivity [1]. In this process the 
released energy Q (equal to the transition energy Qec minus the binding energy of the captured 
electron in the daughter atom) is shared between an emitted electron neutrino and the final atom. 
With a small probability, of the order of 1CT 4 with respect to the normal EC decay, a photon can 
be emitted together with the neutrino. In such radiative electron capture decay (REC) the energy is 
shared statistically between three bodies, thus the energy spectrum of these photons is continuous. 

The theoretical description of REC was initiated by Morrison and Schiff Q2Q who neglected the 
Coulomb field of the nucleus and took into account just IS electrons and only in non-relativistic 
approximation. An advanced and more accurate theory of radiative electron capture was developed 
by Glauber and Martin in Refs. |3i|4|]. They included exactly the Coulomb field in the propagation 
of the electron, took into account relativistic effects as well as screening and considered captures 
from higher shells. However, they limited themselves to the allowed decays where the nuclear 
spin changes by A J = 0, 1 with no parity change. Predictions of this model were tested in a large 
number of experiments, and satisfactory agreement was found yj]. 

A more general theory of REC, extended to any order of forbiddenness (i.e. for arbitrary change 
of nuclear angular momentum and parity), was developed by Zon and Rapoport in Ref . [5] and Zon 
in Ref. Q6Q . For verification of their results, first forbidden unique transitions (A J = 2, 7Tj iif = 
— 1) are of special importance because of a cancellation of the nuclear matrix elements in the ratio 
of radiative to the nonradiative capture rate. Measurements of radiation accompanying the IS EC 
decay in the case of 41 Ca fcfi\, which belongs to this category, revealed a serious disagreement with 
results of Refs. |5i |6fl. The shape of the photon spectrum differed from the prediction, and the 
total probability of the REC process per ordinary nonradiative decay was found to be larger by a 
factor of 6 than the predicted one. To resolve this discrepancy a possibility of photon emission by 
the nucleus, in so called detour transitions, was examined by Kalinowski et al. Q8j, |9|] following 
ideas developed by Ford and Martin [10]. According to Refs. Jsl Hljlj], the nuclear contribution 
to the REC process accompanying forbidden transitions can be substantial. In particular, for the 
case of 41 Ca it was claimed that the detour transitions fully account for the missing intensity 
established by the experiment. However, for another first forbidden unique transition — the IS 
EC decay of 204 T1 — a different situation was encountered. The measured intensity of the REC 
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spectrum ill ill was found to be smaller by a factor of 4 than the value predicted by the model of 



Zon leaving no room for the nuclear contributionjLn contradiction with results of Ref. 

Although we agree in general with Refs. J^l LL(]] that nuclear radiation takes place, we point 
out that the separation of radiation emitted by the electron and the nucleus, respectively, is not 
physical because it depends on the particular gauge used in the description of the electromagnetic 
field. We argue, that although physical results do not depend on the selected gauge, the so called 
length gauge is preferred for the actual calculations. First, it suppresses the nuclear contribution, 
and second, it makes possible important simplifications in the calculations. We note that some 
formulae of Zon in |6|] for A J = 2 transition are divergent for the point nucleus. It means that 
approximations which lead to these formulae may not be correct. For example, an assumption that 
the radius of the region where the photon emission occurs is much larger than the dimension of 
the nucleus where the capture takes place, is not valid in the Coulomb (velocity) gauge used by 
the authors of Refs. flffE)]. 

The failure of the theory of Zon and Rapoport to describe the experimental data for forbidden 
EC transitions, and the controversy concerning the nuclear contribution to the REC process, moti- 
vated us to address these questions theoretically in an independent and different way. In this work 
we recalculate the radiative electron capture process in the length gauge. We restrict ourselves 
to the IS electron (K-capture) only, but the extension to other states is straightforward. In the 
following section we describe our calculations in detail. The results, in terms of a dimensionless 
shape factor, for the case of the first forbidden unique transitions are presented in Section III and 
are compared to the experimental data for 41 Ca and 204 T1. Comparison to previous calculations is 
made in Section IV followed by a summary in Section V. 

II. THEORY 

A. Preliminaries 

In the following we use natural units h = c = 1, e 2 = An a and set the electron mass m — 1. 
The Hamiltonian for the EC-decay is 

H EC = ^ [4 7m (1 - A 7 5 ) ^] [Vv (1 - 7 5 ) V>e] , (1) 

where ip n , ?p p , ?p u and ip e denote neutron, proton, neutrino and electron bispinors respectively, 
with Fermi constant G F = 1.166 39(1) 1(T 5 GeV 2 and A = 1.26992(69). We use convention of 



Bjorken and Drell H12H for Dirac matrices with 7 & = 75 = % 7° 7 7 7 , namely 



'-o-o-o 

and a = 7 7. 

It is convenient to use states with a definite angular momentum. For this reason we introduce 
spin spherical harmonics x™ II 1-411 



re+l/2-m T^m— 1/2 

- I l K+l/2+m m +l/2 I ' (J) 

2k+1 ■ r |K+l/2|-l/2 



where 



j + 1/2 forj = /-l/2 
-(j + 1/2) forj = / + l/2 



(4) 



They have the following properties 



B ■ rC = -X™„ , (5) 
a-Lx™ = + (6) 

which are used to solve the Dirac equation. We consider at first the left handed neutrino, 

i(l-7 5 )Vv = ^- (V) 
Its wave function ty v with the angular momentum j and the energy q is 

= ( j ' (8) 

" = 9 bi-i/ 2 (gr) X m - + i/2)(^ 0) +ijj+i/2(qr) xf+ifitf, 0)] , (9) 

with ji(x) being the spherical Bessel function and v solution of a-pv — — q v. The above neutrino 
wave function is normalized with respect to the energy, namely 

d 3 rip+(f)^ u (r) =27i5(q~q') . (10) 

The electron wave function is the solution of the Dirac equation in the Coulomb field —Za/r 



II 1 311 . It is of the form 

iF{r) xH^M) 



, G(r)x«(M) , 

M^= , (ID 
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with the energy E = E n)K 



E. 



n,n 



Z 2 a 2 



1 + 



> + I + 7) 2 



-1/2 



where 7 = — 1 + V k 2 — Z 2 a 2 . The explicit form of radial wave functions G and F is 1131] 

G n>K (r) = C„,^2Arre- Ar v ^T^:^S(r)+^)(r)], 
F n>K (r) = C nt 42Xrye- Xr ^^[gZ(r)-g^(r)}, 
gM(r) = VZa\-i-KL^\2\r), 



,(2), 



= + 2 + 2 7) VZ a A-i - « L^ 2 7) (2 A r) 

G n ,K 



2A 4 n! 



Zar(n + 3 + 2 7 ) 



(12) 

(13) 
(14) 
(15) 
(16) 

(17) 



where L„ are Laguerre polynomials 



fc=0 



k\ (n- k)\T(a + k + 1) 



-z 



(18) 



with L-i = and A = \/l — F 2 _ In particular, the energy and the wave function of the ground 
state n — I — 0, k — — 1, and j — 1/2 are given by 



#o,-i = Vl-^2 a 2 =£ 



/4Z 3 a 3 (2 + 7) 

r(3 + 2 7 ) 

'4Z 3 a 3 (-7) 



(2Zar) 7 e 



7 „—Z a r 



0,-1 



(2Zar) 7 e- Zar Ef. 



(19) 
(20) 

(21) 



r(3 + 2 7 ) 

For the calculation of radiative transition rates one needs the Dirac-Coulomb Green's function G D . 



Its explicit form is fl 1 3M 

G D (f,f',E) = (f 



H-E' 



n = J2 G ^m(r,r',E) 



Gl m {r,r',E) = 9(r' - r) $<J?) ® ^(r') 1 + e(r - r') V£ m (0 ® ^(f) 1 
where 0(x) is a Heaviside step function, and both and ^ < are of the form (fTTI) with 



G<(r) = (2Arre- Ar v / r+#(/ 2 + A) 



F<(r) = (2Arre- Ar v / T^#(/ 2 -/i), 
G>(r) = (2\rye-* r VTTE(U + f,) 2 -^- 

r(c) 

F>(r) = (2Arre- A ^v / T^(/4-/3)^|^ 



(22) 
(23) 

(24) 
(25) 
(26) 

(27) 



where 



fi = (ZaX 1 — k) i-Fi(a, c, 2 A r) , 

/ 2 = a 1 F 1 (a + l,c, 2Ar) , 

/ 3 = f/(a,c,2Ar), 

/ 4 = (ZaA- 1 + /t)[/(a + l,c,2Ar) 



(28) 
(29) 
(30) 
(3D 



and a = l + 7 — E Z a X 1 , c = 3 + 2 7, while 1-F1 and Z7 are confluent hypergeometric functions 
regular at the origin and at the infinity, respectively. 



To describe photon wave function we introduce vector spherical harmonics II 1411 . 



Yfdo, 4>) = zZ <t>) ^ ( L ™; 1 <?|£ i; JM ) 



(32) 



m q 



where e\ = — 1/ \[2 (e x + ie y ), e = e z , and e_i = l/v2 (e x — ie y ). The solutions of the Maxwell 
equations with definite angular momentum and parity are represented by both, the magnetic photon 



,0(M) 



'2*jj(*r)y#(r 



(33) 
(34) 



and the electric photon in the length gauge fl 1 3M 



21^l 2 -^J J+1 (kr)Yf J+1 (f), 



A°£\r) = -iV2k\/^jj(kr)Yf(r) 



while the electric photon in the Coulomb (velocity) gauge is 

V2k 



(35) 
(36) 



^-j-py jj+i(k r) Yj§ +1 (r) -J^Ljj^kr)?"^ 



0. 



(37) 
(38) 



These solutions are normalized with respect to energy, so that 



2 k <Pr \Ajm(j)* A' JM (r) - A° JM (r)* A%(r)] = 2n5(k- k') 



(39) 



The use of the length gauge is essential for performing several simplifications in the calculations 
of transition rates, as it is explained in the next sections. 



B. Electron capture rate 



The electron capture rate W is equal to the square of the the matrix element M, summed over 
final states and averaged over initial states 



W 



2 J, 



— T m 



(40) 



M e Mi M„ M/ 



where 



Gi 



M = (f\H EC \t) = J d 3 r-^ llx (1 - A 7 5 ) Vi] 7 M (1 - 7 5 ) • (41) 

Although we use single nucleon matrix elements, results can easily be transformed for the nuclear 
matrix elements, by assuming that ip n and tp p are field operators, and instead of ipf 7^ (1 — A 7 5 ) ipi 
one considers (f\ipn1^ (1 — A 7 s ) -0_p|z). The tensor decomposition of the matrix element in Eq. 
(I4TT) leads to 

J JLSM 



where 



T 



T 



M 
JLO 

M 
JLl 



i L Sjl Y L M , 



(43) 
(44) 



and (.|.) denotes the integral over angular coordinates. Each state f,i,ip v and 4> e has definite 
angular momentum J, M numbers, so one can use reduced matrix element 



(j,™\ T Z\j ,m) = (-1) 3 ' 



J k f 



-m q m 



(45) 



and orthogonality properties of 3j symbol Ill4h to obtain simple formula for the electron capture 
rate 

■»2 



W 



G l F 1 1^1 

~Y 2 j e + 1 2 Ji + 1 ^ 2 j + 1 

/ r 2 dr ^ (J/||T JL5 (1 - A 7 5 )|| J t ) (J V \\T JLS (1 - 7 5 )|| J e ) 



(46) 
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The reduced matrix elements of spherical harmonics are given by 



(*/IMI«i> = y^^C^/,^), (47) 

Cfaf,*) = (-l)^ +1 / 2 V / (2j / + l)(2j l + l) I 3t 1 3f )u(l f ,k,l), (48) 

v y 1/2 -1/2 y 

n(l/, Ig, Z) = ^[l + (-l) l f +k+l ], (49) 

and spin spherical harmonics by 



(Kf\\Y JL ■ a\\iu) = Y 4^ 1 S JL {K, f ,K,i) , (50) 

Sj,j +1 (Kf,Ki) = ( 1 + ~j+t) Cj ^~ Kf ' ' (51) 

Sj,j(Kf,Ki) = -^====L= Cj(K f ,Ki) , (52) 

*) = ^TTT (-1 + C7 J (-« / , *) . (53) 

With the use of the above formulae for reduced matrix elements, the capture rate W of the IS 
electron by the nucleus is 

W = J2 W J> (54) 

j 

^° = |x4^4^ |(J/ll[T000 ^ ( ^" T011 ^ )](1 " A75)l|Jt)|2 ' (55) 

(56) 

W2 = l^|fT4^l <J/l|r [ T 2ii^(i?) + (T 220V ^ + T 221% /375)^(i?)](l-A 75 )||J i )| 2 , 

(57) 

where Q is the energy released in the decay and R is the nuclear radius. 
C. Radiative electron capture rate 

The probability amplitude for the electron capture with the simultaneous photon emission is 

M R = J d 3 r^= 7 M (1 - A 7 5 ) A] [Vv Y (1 - 7 5 ) €\ , (58) 
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where 

n 

We use the latter form to perform tensor decomposition of M R and obtain 
Mr = —7=- I r",ir > ( — 1 } J + A] + l (,/, . M r \'.rj[ s ( I - A V* ) I ■/, M ) > (60) 



v/2 



fr-dr (-1) J+M+1 (^^/I^5(1-A7 5 )I^M0 £ 

JLSMM' 



(J„, M^TJ^ (1 - 7 5 )|n, ^, ikQ En + k _ g (n, J' e , M>\A°j a m a ~ « ■ A/.mJ J b , M.) 
The rate for the radiative electron capture H^j is 

W fl = / rffc W R (k) , (61) 







where 

W R (k) = V — \M R \ 2 , (62) 

v ; 2 J e + 1 2 Ji + 1 ^ 2vr 1 ' 



MeMiMuMfM A 



with the sum over final state and the average over initial states. The W R (k) can be expressed in 
terms of the reduced matrix elements and the summation over magnetic states can be carried out 
with the aid of the orthogonality of the 3j symbols. The results is 



w R( k) = o- , 2 1 2 1 ^ s ^ie^z 

x / ^^(J/HWi-At 5 )!!^) (J,||T Ji5 (i- 7 5 )||^,n) 

J 7-0 



(63) 



We can use now the explicit form of the Dirac-Coulomb Green's function in Eq. (I22H31I) to replace 
the sum over intermediate electron states n in Eq. (1631 by G D . Since the electron capture takes 
place within the nucleus, and the photon radiation in a region of the electron wave function which 
is several orders of magnitude larger, we apply an identity 0(r' — r) = 1 — 0(r — r') and neglect 
6(r — r') in Eq. (1231) completely, so the Green's function becomes 

G D (r, f E) « ^< m (f ) ® ^> m (f 0+ • (64) 
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In other words, this approximation is allowed, because the integral with 0(r — r') gives contribu- 
tion, which is higher order in the small parameter £ = Q R. After this assumption one obtains 

x / r 2 dr^ (J/HW1 -A 7 5 )||Ji) WIWI-t 5 )!!^^)! 2 

42X^ |(J: '^ >ll ^~ C5 ' ljJ|Je)|2 }- (65) 

The approximation in Eq. (l64l) would not be valid in the Coulomb gauge, as the integral with 
6(r — r') is of the same order in the parameter £, what we discuss in more details in Sec. IV. 
We consider transitions with A J = 0, 1, 2 

W R (k) = W R0 (k) + W R1 (k) + W R2 (k) , (66) 

since no electron radiative capture has been observed for higher multipolarities. Due to the as- 
sumption in Eq. (l64l) . transition rate for each value of J can be decomposed into a product of 
a term corresponding to the photon emission and a term corresponding to the nuclear transition, 
namely 

W m {k) = W M1 (k, S 1/2 - S> /2 ) W (Q - k, S< 2 - u 1/2 ) 

+W E1 (k, S 1/2 -> P> /2 ) W (Q - k, P< /2 -> i/i /2 ) , (67) 

WOu(fc) = 5 1/2 - S> /2 ) Wi(Q - fc, 5< 2 - z/ 1/2 ) 

+Wk(fc, S 1/2 -> P> /2 ) W X {Q - fc, P< 2 -> i/ 1/2 ) , (68) 

iy ra (fc) = W M1 (fc, 5 1/2 - S> 2 ) W 2 (Q - fc, S< 2 - u 3/2 ) 

+W E1 (k, S 1/2 -> P> 2 ) W 2 (Q - fc, P< /2 -> ia, /2 ) (69) 
+ [W E1 (A;, 5 1/2 - P 3 > /2 ) + ^ M2 (A:, 5 1/2 - P 3 > 2 )] W 2 (Q - fe, P< 2 - z/ 1/2 ) 
+ [^mi(A;, S 1/2 - P> /2 ) + W^ 2 (fc, S 1/2 - D> /2 )] W 2 (Q - k, D< 2 - z/ 1/2 ) . 

The explicit formulae for the photon emission rate 

WW = 2j^Y I ^ W A °i A ~ « ■ ^ 1 1 4.) I = - fc Kem (70) 
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are obtained using photon wave function from Eq. (T33H36I) . electron wave function from Eqs. 
(12012 II) . and the Dirac-Coulomb Green's function from Eq. (l64l) 

2 



7^a/i(Si/2 - > S 



1/2 ) 



r 2 d rjl (kr) [G> x {8 - k,r) F{r) + F> x {£ - k,r)Q(r)] 
: dr j x {kr) [G>(£ - k, r) F(r) + F>{£ - k, r) Q{ 



T^Ml (S1/2 —> -D3/2) 



n M2 (S 1/2 -> P>) = 3 / r 2 drj 2 (kr) [G> 2 (£ - k,r) F{r) + F> 2 {£ - k,r)Q{i 



fcEl(Sl/2 — > Px/ 2 ) 



fcEl(Sl/2 -> P3/2) 



7lE2(Sl/2 — > -D3/2) 



y r 2 rfr|j!(A;r) [Gf (£ - k, r) Q{r) + F?{£ - k, r) F{i 
+2j 2 (kr)G>(£ -k,r)F(r)}\ 2 , 

y r 2 rfr |2ji(fcr) [G> 2 (£ - k,r) Q(r) + F> 2 (£ - k,r) F(i 
+j 2 {k r) [G> 2 (£ - k, r) F{r) - 3 F> 2 (£ - k, r) Q{i 



J r 2 dr{j 2 (kr) [G 2 (£ - k, r) Q{r) + F 2 (£ - k, r) T{i 



+2j 3 (kr)G>(£-k,r)F(r) 



(71) 



(72) 



(73) 



(74) 



(75) 



(76) 



where E = Eq-i, T{r) = Fo ) _i(r), and G(r) = G*o,-i( r )- The explicit formulae for the electron 
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capture rate are 



W (S< /2 - u 1/2 ) = 2 -^f^ ^\(J f \\[T 000 G< 1 (£-k,R) 

-Ton F<,(£ -k,R)](l-X l5 )\\J t )\ 2 , (77) 

W (P< 2 ^u 1/2 ) = 2 °\ f'^ ^ I (41 1 [Tooo F<(£ - fc, R) 

+T 011 Gf(£ - k,R)] (1 - A7 5 )|| Jj)| 2 , (78) 

Wl (S< 2 ^u 1/2 ) = 2 G \^~^ f ^ | (Jf\\ [T 10 i gf t (g - fc, R) 

+ {T 1W /V3 + T 1U v^73) F< (£ - k,R)] (1 - A T 5 )||^>| 2 , (79) 

W l (P< 2 ^ u l/2 ) = 2 -^ff^^\(Jf\\[TioiF<(S-k,R) 



(T 110 /V3 + T m v^73) G<(£ - fc, i?)] (1 - A 7 5 )|| Ji)| , (80) 



^ 2 (S< 2 ^ 3/2 ) = ? G \ { f i + ^ IW G<i(g ~ fc, A) 

+ (T 220 v 7 ^ + T 221 v^75) F< X {E - k,R)] (1 - A7 5 )||J;)| 2 , 

W 2 (P< /2 - i* /2 ) = I ^ ^ |(J/lir [T 211 F<{£ - k, R) 

- (T 220 + T 221 v^75) Gf {£ - k, R)] ( 1 - A 7 5 ) 1 1 J*) f , 

a 2 (o - k) 2 1 

W 2 {P< 2 ^v l/2 ) = {2 F j W +1) ^ ^\(Jf\\r [T 211 G< 2 (£-k,R) 

VV5 + T 221 v^S) F< 2 (£ - fc, ( 1 - A 7 5 ) 1 1 Ji) \ , 

W 2 (D< /2 ^ u l/2 ) = fj [Q + ~^ 2 ^\{Jf\\r [T 211 F<(S-k,R) 

- (T 220 V2j5 + T 221 G< {£ - fc, R)] (1 - A 7 5 ) 1 1 J,) | 2 . 



(81) 



(82) 



(83) 



(84) 
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III. RESULTS 

Following the generally adopted convention Jsl, O, Isl, we calculate the ratio of radiative to 
nonradiative capture rate. Then, the resulting photon spectrum accompanying the IS capture is 
given by 

W R (k) __ a ( k x 2 



1- kTZ, (85) 
W nm z \ Q J 

where details of the shape are expressed in terms of the dimensionless shape factor 1Z. 

First, we consider the transitions where A J = 0,1. When the parity of nuclear states does 
not change, Hi iif = 1, this transition is allowed. Then, two tensor operators T 000 and T 101 (— A7 5 ) 
dominate and other operators can be neglected. In the probability ratio the nuclear matrix elements 
cancel out and the result greatly simplifies: 



H - Ha 1 gT^W — ^ K M i{Sx/2 — > Sy 2 ) + 1 ^ g(^ R y — ~ ^Ei(Si/2 -> Pi/ 2 ) ■ (86) 



It is worth to note, that in the nonrelativistic limit (Z a — > 0), the shape factor 1Z a is equal to 1, 
what we verified by explicit calculations. If the parity does change, 7Tj Hf = —1, the transition is 
first forbidden non unique. Here, all nuclear operators contribute to the transition matrix elements 
and the exact expression for 1Z in general is more complicated. However, one notices that the ratio 



F<(E,R) [T^E VK 2 -Z 2 a 2 + K-Za^/]±§ (2k\ kM G< k {E,R) 

(87) 



G<(E,R) \l + E^ K2 _p Q2 _ K + Za f±E \ZaJ F< K (E,R) 



l+E 

approximately does not depend on energy E as long as Z a 1 . Therefore, nuclear matrix 
elements cancel out and the formula in Eq. (f8~6l) holds for non unique transitions, but only when this 
approximation is valid. It is important to stress that numerically Eq. (|8~6l) yields the same results 
as those obtained in Refs. J3], 0, Q]. Although our results have apparently different and simpler 
form, they are equivalent to those obtained previously by Glauber, Martin, Zon and Rapoport for 
the allowed and first forbidden non unique decays. 

The situation is different, however, in the case of the first forbidden unique transition, where 
A J = 2 and 71^717 = —1. The operator T 2 n(— A7 5 ) dominates and other nuclear operators are 
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neglected. The resulting shape factor is 

A £;\ 2 _^ k 2 

rClu 



Q 



1* 7^(i) + 4_ 7^(2) 



n (2) 



9 G%{£-k,Rf 



2k 2 R 2 G(R) 2 
4 9 F<(£-k,R) 2 



[7^ei(5'i/2 - 
[7^- Ail (51/2 



-^3/2) 



2k 2 R 2 g(R) 2 
While the nonrelativistic limit of IZ^ is 1, for the IZ^ function we obtain: 



^3/2 



U Zj2 



lim Tl {2) 
z->o 



1 2 

1 + 



(88) 
(89) 

(90) 
(91) 
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what we have verified numerically. Although the result of Zon [6] has similar form to Eq. (f88~T) . the 
Coulomb-free limit of the function corresponding to IZ^ is equal to 1, which differs significantly 
from our result. 

Experimentally, the REC process in the first forbidden unique transition was studied so far for 
two nuclei: 41 Ca [7] and 204 T1 JuJ]. The data obtained for 41 Ca (Q EC = 421.3 keV, T 1/2 = 10 5 
years) are shown in Fig. 1 together with theoretical predictions. While the model of Zon and 
Rapoport [3, 6] underestimates the intensity of the REC spectrum by a factor of about 6, the results 
of our calculations are in excellent agreement with the experiment. The contribution of relativistic 
and Coulomb effects is evident. The Coulomb-free limit of our result would overestimate the 
intensity approximately by a factor of 2. The discrepancy between the Zon model and the 41 Ca 
data was used as an argument in favor of nuclear contribution to the REC process [^l 9]. In fact, if 
both, the nuclear and the electron radiation were calculated correctly, the result should agree with 
Eq. (|8~8~T) . what is explained in the next section. 

The case of 204 T1 (Qec — 346.5 keV, Tu 2 = 3.78 years) is illustrated in Fig. 2. Again, our 
results perfectly agree with experimental data, while the results of Zon and Rapoport overestimate 
the spectrum by a factor of about 4 111 ill . As expected, the influence of Coulomb field is much 
stronger in case of thallium (Z — 81) than in case of calcium (Z = 20). The nonrelativistic limit 
would overestimate the intensity of the REC spectrum by a factor of about 10. 



IV. COMPARISON TO PREVIOUS CALCULATIONS 

We have verified numerically that Eq. ([86]) agrees with former results obtained in Ref. [|3i|4|,|5, 
6|] for the allowed and first forbidden non unique EC transitions. However, in contrast to previous 
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FIG. 1: The shape factor for the IS component of the first forbidden unique REC in 41 Ca. The points with 
error bars represent the experimental values form Ref.ITj]. The solid line shows the result of this work as 
given by Eqs. 86-88, while the nonrelativistic limit is shown by the dashed line. The dotted line represents 
the prediction of the Zon and Rapoport la,|6D- The value of Q is 416.4 keV. 



works [51 la |9fl, our results correctly describe the known experimental data for the first forbidden 
unique decays. The difference, in our opinion, is due to incorrect approximations adopted in 
the Refs.o, |6J]. If we use the Coulomb gauge for the emitted photon, Eqs. (1371381) . then the 
approximation of Eq. (|64l) for the Coulomb Green's function is not valid, because the integral 
with the neglected remainder, which is proportional to 0(r — r') is as important as the integral 
with the approximated form in Eq. (l64l) . As an consequence of this approximation in the Coulomb 
gauge, 1Zei{S\/2 — > -P3/2) in Eq. (1751) would contain jo(kr), and the integral would diverge for 
small values of r, as the integral in the function B 2 \ of Ref. [6] does. With the length gauge, Eqs. 
(1351361) . where El photon wave function contains a combination of ji (k r) and j 2 (k r) functions, 
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FIG. 2: The shape factor for the IS component of the first forbidden unique REC in Tl. The points with 
error bars represent the experimental values form Ref.[ll]. The meaning of theoretical lines is the same as 
in Fig.l. The value of Q is 264.0 keV. 

the r-integral in Eq. (1751 ) is finite, as it should. Moreover, the probability amplitude for the 
radiation from the nucleus is strongly suppressed, because the photon wave function is much 
smaller within the nucleus. However, we can not exclude completely the nuclear contribution if 
there are close-lying excited states in the daughter or in the parent nucleus. In such a case, the 
nuclear radiation should be obtained using the same gauge. Since the photon wavelength is much 
larger than the size of the nucleus, the dipole approximation is allowed. Then, the coupling of 
nucleus to the electromagnetic field in the length gauge takes the form —d • E — Jl ■ B, and this 
form should be used for the calculation of the photon radiation from the nucleus. 
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V. SUMMARY 



We have recalculated the rate of radiative orbital electron capture by the atomic nucleus. We 
applied the length gauge for the radiated photon, because it suppresses the nuclear radiation and 
substantially simplifies the calculations. By the use of a convenient form of the Dirac-Coulomb 



Green function from Ref. 111311 . the results for the capture of the IS electron can easily be gen- 
eralized for the arbitrary state of the captured electron. The results obtained are found to be in 
good agreement with experimental data for allowed, first forbidden and first forbidden unique 
transitions, in particular with the results for 41 Ca and 204 T1. 

There is another case of the first forbidden unique EC transition, the decay of 81 Kr, which can 
be used to verify predictions presented in this paper. An experiment with a sample of 81 Kr is being 
undertaken, and its outcome will be published soon. 
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